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ABSTRACT 

It has long been known that the energy in velocity and magnetic field fluctuations in the solar wind is not 
in equipartition. In this paper, we present an analysis of 5 years of Wind data at 1 AU to investigate the 
reason for this. The residual energy (difference between energy in velocity and magnetic field fluctuations) 
was calculated using both the standard magnetohydrodynamic (MHD) normalization for the magnetic field and 
a kinetic version, which includes temperature anisotropics and drifts between particle species. It was found 
that with the kinetic normalization, the fluctuations are closer to equipartition, with a mean normalized residual 
energy of CTi = -0.19 and mean Alfven ratio of ta = 0.71. The spectrum of residual energy, in the kinetic 
normalization, was found to be steeper than both the velocity and magnetic field spectra, consistent with some 
recent MHD turbulence predictions and numerical simulations, having a spectral index close to -1.9. The local 
properties of residual energy and cross helicity were also investigated, showing that globally balanced intervals 
with small residual energy contain local patches of larger imbalance and larger residual energy at all scales, as 
expected for non-linear turbulent interactions. 
Subject headings: magnetic fields — MHD — plasmas — solar wind — turbulence 



1. INTRODUCTION 

The solar wind displays many properties that are consis 
tent with our current id eas of plasma turbulence (see re 
cent reviews bylBruno & Carbone 2005; Horburv et al. 200^ 
Petrosvan et al I 120 lOt iMatthaeus & VeUi, ,2011.: .Carbont . 
2012h . These include power sp ectra of field fl uctua- 



I 



tions that have power law forms (Colemad 119681) . non 
Gaussian probability density functions (PDFs) of field incre 
ments (Marsch & Tu 1994) and various forms of anisotropy 
JHorburv et al...201 1 ). One observational feature that is less 
well understood, however, is that the energy in velocity and 
magnetic field fluctuations is not generally in equipartition. 
In this paper, we present new measurements to investigate the 
reason for this. 

The power spectra of magnetic field and velocity fluc- 
tuations in the solar wind are generally interpreted in 
terms of magnetohydrodynamic (MHD) turbulenc e the- 
ory (M roshnikov 1963; Kraichnan 1965; Goldreich & Sridhai 
1995nGaltier et al.ii2000E iBoldyrevii2006t iSchekochihin et al. 



for scales larger than the ion gyroradius. Since such 
theories involve the non-linear interaction of Alfvenic wave 
packets, the energy in the velocity and magnetic field fluc- 
tuations has been assumed, either ex plicitly or impl icitly, to 
be in equipartition, as it is for pure ' Alfvenl (Il942h waves. 
However, it has long been known that the energy in mag- 
netic field fluctuations dominates the energy in velocity fluc- 
tuations (e.g.. Belcher & Davis 1971; Matthaeus & Goldstein 
1982; Bruno etal. 1985; Roberts et al. 1987; Tu et al. 1989[ 
Marsch &Tu 1 990; .Gr api^in et al. 1991; Goldstein et aTl 
1995; Bavassano etal. [1998; Podesta et al. 2007; Salem et al. 



iPerri & Balogh 



120 lOt iChen et al.i i2011bt iBorovskvi 



Furthermore, the spectral index of the magnetic field and 
velocity fluctuation s has been measured to be different. 
iGrappin etal.l(ll991h showed that the velocity spectra are sys- 
tematically shallower than the magnetic field spectra, which 
has been confirmed with more recent measurements show- 
ing an average velocity spectral index close to -3/2 and mag- 



netic field spectral index close to -5 /3 (Mangenev et al. 2001 ; 
[Podesta et al.' 120071; 'Tessein et all [2009 ; Salem et al. 2003 
Podesta & Borovskv 2016; Ch en etal.1 1201 la; Boldyrev et aL| 
[20 Ik lBorovskyll2012l) . This difference has also been con- 
firmed with electric field measurements (Chen et al. 201 la). 

An excess of magnetic energy is frequently, but not 
always, seen in numerical simulations of MHD turbulence, 
both forced and decaying, weak and strong, balanced 
and imbalanced and with and without a mean mag 
netic fi eld (e.g., lOughtonetal.1 [l994; Biskamp & Mulle 
19991: iBeresnvak & Lazariah 2006; Muller & Grappi" 
20051 iBTgot et al. 2 008; Boldvrey & Perez 2009L 
Mininni & Pouquet 2009r iWang et all [20111: IChenetall 
i2011bi) . The magnetic field spectrum is also often 

steepe r than the veloc ity spe ctrum in MHD simula- 

tions (Muller & Gra ppinI [20051: [Mininni & Pouquej [2009[: 
Boldyrev etal. 2011)' 

Several possibilities have been proposed to explain these 
differences between the magnetic field and velocity fluc- 
tuations in the solar wind. Early explanations involved 
the neglect of non-MHD corrections to the Alfven speed 
from temperature anisotropics, whi ch would affect the nor- 
m alization of th e magnetic field dBelcher & DavisI [19711: 
Matthaeus & Goldstein 19 82). Other aufliors (e.g., [HoUwegl 
1987; Tu et al. 1989; Robe rts et all 19 90), however, suggested 
that the typical measured anisotropies were not suffi cient to 
explain the difference and [Bavassano & Brund ([20001) showed 
that when incorporating the anisotropies, the magnetic fluctu- 
ations, while sometimes smaller, still dominate. 

An alternative explanation is that the difference between 
the energy in velocity and mag netic fi eld fluctuations, called 
residual energy (Pouquet et al .|[ 19761), is an inherent feature 
of the turbulence itself. Grapp in et al.l ([19831) showed numer- 
ically and analytically that turbulence in an isotropic closure 
theory has residual energy with a wavenumber spectrum k~^, 
independent of the level of imbalance between fluxes of op- 
po sitely directed Alfven wave packe ts. This was expanded on 
by [Muller & GrappinI ([2004 [2005[) , who showed, using the 
same closure theory and MHD simulations, that the residual 
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energy spectrum would be k ''l^ for an energy spectrum of 
k~^l^ and k~^ for an energy spectrum of k~^l^. They described 
the presence of a small residual energy as a competition be- 
tween magnetic field amplification from a local dy namo effect 
and e quipartition from the nonlocal Alfven effect (iKraichnanI 
fT%5h. 

iBoldvrev & PerezI (120091) proposed that weak turbulence 
naturally generates a condensate of residual energy in small 
parallel wavenumber (A:|| = 0) modes (which do not not 
have to satisfy the equipartition since they are not waves) 
due to mirror-invari ance breaking in imbalanced turbul ence 
(see also H chekochi hjn et al.ll2012h . IWang et alJ (1201 Ih and 
iBoldvrevet al., (,20 12b *) showed analytically and numerically 
that this residual energy is statistically negative, i.e., the mag- 
netic energy dominates, and has a perpendicular spectrum 
of k~^ . Using an anisotropic closure theory, Boldvrev et al.l 
( I2012al) predicted strong turbulence to have a k~^ resid- 
ual energy spectrum, independent of the particular strong 
turbulence model. Simulations of forced, strong reduced 
MHD turbulence show a scaling of fc~ '^ for both balance d 
and moderately imbalanced turbulence (Soldvrev et al.l20T11) . 
iGogoberidze et al.l (l2012al) also used an anisotropic closure 
theory, assuming strong, balanced, axisymmetric turbulence, 
finding the residual energy to be negative but with a k~^l^ 
spectrum. 

Further suggestions for the residual energy, that are not 
necessarily contradictory to the above explanations, in- 
clude the presence of a separate component of convect- 
ing magneti c structures, i n wh i ch th e re are no velocity 
fluctuations (Tu & Marsch' 1991, 1991 iBruno et all |2007|) . 
iMininni & Pouquet (2009) used simulations to show that the 
different scaling of the velocity and magnetic fields could be 
due to their different amounts of intermittency, the magnetic 
field being more intermittent due to its ability to form thin cur- 
rent sheet s, and therefore having a steeper power spectrum. 
Similarly, iLi et al. (2011) suggested that the magnetic field 
takes a -3/2 spectrum to match the velocity in the absence of 
current sheets, although whether such discontinuiti es are part 
of the turbulence or not remains an open questio n (iBorovskvl 
I2008t iGreco et al.ll2008t IZhdankin et al."2012b^. Compres- 
sive fluctuations have also b een sug gested to be a cause of 
the residual energy, although Hollwegj (119871) noted that their 
presence would lower, rather than increase, the magnetic en- 
ergy. 

While the amount of residual energy has been frequently 
measured in the solar wind, its scaling properties have been 
rarely documented. To our knowledge, the only measure- 
ments of residual energy scaling were done by Tu et al. 
(|l989), who found a slope near to -5/3 for an interval of fast 
wind measured by Helios at relatively low frequencies (below 
6 X lO""* Hz) and bv lPodesta & Borovskvl (120101) who quoted 
a value of -1.75 but gave no further details. In this paper, 
we present a large survey of residual energy scaling measure- 
ments, compare these to the above theoretical predictions and 
make new measurements of the local properties of the residual 
energy. 

2. DATA SET 

To make reliable scaling measurements of solar wind tur- 
bulence, it is necessary to use a large amount of data. While 
it would be better to use a single long interval in which 
to measure the spectral indices, in practice this is not pos- 
sible in the solar wind, so an average of the spectral in- 



dices from many shorter intervals is typically used. His- 
tograms of inertial range spectral index measurements made 
in small intervals, typically a few hours long, show that 
they are approximate ly normally distributed with values rang 



ing from -2 to -1 (Smith et al. 2006; Vasquez et al. 2007b; 
ISssein et al. 2009; Chen et al. 2011a; Boldvrev e t al. 201 1; 
Borovskv 2012). Simulati ons also show a simil ar spread of 
values between snapshots (IBoldvrev et al.l |201 Ih suggesting 
this to be an inherent feature of the turbulence. Although 
the mean of this distribution does not have to be the same 
as that calculate d from a single l ong interval, observat ional 
( iBorovskvl2012l) and simulational (IBoldvrevet al .1120 lib data 
suggests that they are. 

For the current analysis, 5 years of data from the Wind 
spacecraft (Acufiaet al.lll995l) . from June 2004 to April 2009, 
were used when the spacecraft was continuously at the first 
Lagrange point (LI), ?s230 Re upstream of Earth. This period 
covers the declining phase of solar cycle 23 and the start of 
solar cycle 24. Magnetic field data from MFI (Lepping et 
19951) . onboard ion moments from 3DP /PESA-L dLin et 



I995I) and ion parameters from SWE/FC (lOgilvie etal.ll 19951) 



were used in the analysis. 

The data were split into 6 hour intervals in which various 
parameters, including the spectral indices, were calculated. 
This interval length is long enough to enable a reliable spec- 
tral index to be measured over our chosen frequency range 
but short enough so that there is minimal mixing of different 
solar wind types and the variation of the scaling with back- 
ground parameters can be determined. Small data gaps in the 
magnetic field data and proton velocity moments were lin- 
early interpolated over and intervals with more than 5% data 
gaps or containing less than 5.4 hours of data were excluded. 
Several intervals in which 3DP was not in the correct mode 
(during January 2005, November 2006 to January 2007 and 
August to September 2007) were also removed, leaving 5,990 
intervals of 6 hour duration. 

We have chosen not to remove any discontinuities, sudden 
changes of the magnetic field, from the data. It is perenni- 
ally debated whether these discontinuities are generated by 
the turbulence or are convected static structures related to pro- 
cesses a t the Sun (e.g., Burlaga 1969; Neugebauer et al. 1984 
Horburv et al. 2001; Bruno et al. 2001 : Vasquez et al. 2007ai) . 



Borovskvl (12008!) argued that the strong discontinuities are 
due to convected flux tubes and IBorovskvl (120 1 0*) showed that 
they contribute significant power to the fluctuation spectrum. 
ILi et al.l (120111) used a few intervals free of current sheets to 
suggest that the magnetic field spectrum takes a -3/2 scaling 
in their absence, although Borovsky (2012) found the oppo- 
site trend, with the magnetic field spectral index becoming 
shallower for larger numbers of strong discontinuities. On the 
other hand, it has been shown that turbulence simulations can 
reproduce many of the properties of the solar wind disconti- 
nuities, such as the PDFs of waiting times (Greco et al. 20081 
I2OO9I) and field rotations ( Zhdankin et al...2012b..a) . Borovskvl 
(2010) showed that the strong discontinuities alone produce 
a -5/3 spectrum, which could be interpreted as them being 
generated by the turbulence. In this analysis, we have cho- 
sen not to remove any discontinuities from the data, since the 
majority are plausibly part of the turbulence. 

The magnetic field, B, from MFI and the proton velocity, 
V, from 3DP/PESA-L were used to construct three further 
quantities: the magnetic field in velocity units, b, and the two 
'Elsasseii (ll950l) variables, z^ = v±b. Since z^ correspond to 
Alfvenic wave packets traveling parallel and anti-parallel to 



B dKraichnanll 1 965h . the direction o f B was "rectified" based 
on the predicted local lParkerl (11958 1) spiral direction so that z^ 
corresponds to outward Alfvenic propagation away from the 
Sun and z" corresponds to inward propagation towards the 
Sun dBruno et al.l[T98 5l 'Roberts et al.| |1987D . 

The normalization of b is a subtle issue. Usually the nor- 
malization 



B 



(1) 



is used, where p is the mass density (based on the proton den- 
sity or with an estimated fraction of alpha particles). This 
makes the incompressible MHD equations symmetric in z^ 
and we call this the "MHD normalization." For a kinetic 
plasma that is not in thermal equilibrium, such as the solar 
wind, the Alfven wave remains noncompressive but is modi- 
fied, with corrections for temperature anisotropics and relative 
drifts be tween species . One can then normalize the magnetic 
field as (iBarnesll 19791) 



b = 



B 



1 + 



\p^-p\\-y^m,n,{l\ysf 



(2) 



where p is the total mass density, p±^ and p\\ are the total per- 
pendicular and parallel pressures, m^ is the mass of species 
s, lis is the number density of species s and Av, is the drift 
of species s with respect to the centre of mass velocity. We 
call this the "kinetic normalization." In Section [3] we discuss 
results using both normalizations. 

For the MHD normalization, the mean p for each 6 hour in- 
terval was calculated in the standard way, using proton num- 
ber density data from 3DP/PESA-L and assuming that 4% of 
the solar wind ions are alpha particles, the rest being pro- 
tons. For the kinetic normalization, proton and alpha pa- 
rameters derived from SWE measurements (as described in 
[Maruca et al. 2011; Maruca & Kasper 2012) were used to de- 
termine the 6 hour average values of p, p±, p|| and Av.,. Since 
the electrons are fighter and generally more isotropic than the 
ions (due to their higher collisionality and/or kinetic instabil- 
ities (Stverak et al.l 120081) ). we expect them to have a much 
smaller contribution in Equation |2] they have not been in- 
cluded in the present analysis, but will be investigated in fu- 
ture work. It was also assumed that the other higher mass 
ions are sufficiently scarce to contribute to the normalization. 
An interesting property of the kinetic normalization is that it 
occasionally produces an imaginary magnetic field when the 
plasma is unstable to pressure anisotropy or drift instabilities. 
These intervals, along with those where no SWE data was 
available, were excluded from the analysis, leaving 5,396 in- 
tervals of 6 hour duration. 

In each 6 hour interval, the trace power spectra of v, 
b and z^, denoted Ey, E^ and E±, were calculated using 
the multitaper m ethod with time-bandwidth product NW = 4 
(iPercival & Wald en 1993). The spectra of residual energy E^, 
cross helicity E^, total energy Ei, normalized residual energy 
(Tr, normalized cross helicity a^, Alfven ratio ta and Elsasser 
ratio te, were then calculated. 



Ei:=E+-E-, 

Et = Ey + E\,, 

Ey -Eb rp^- 



1 



cri- = 



Ey + Eb ta+I 



(3) 
(4) 
(5) 

(6) 
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Figure 1. Power spectra of magnetic field (Et), velocity (E,) and residual 
energy {Er = E^ -Eb). The light green line represents negative £1 and the 
dark green crosses represent positive Ey. The solid black lines are fits to the 
measured E\, and £v spectra and the difference of these fits and have slopes 
of -1.66, -1.52 and -1.96 respectively. 



ar = 



E+-E- _ rE-l 

E++E- ~ rE + r 

Ey l+CTr 



^^=^=1 



l+CTc 



r =^ = 

'"^ E- 1-cr, 



(7) 
(8) 
(9) 



While the centre of mass velocity should ideally be used to 
calculate these spectra, the alpha, electron and minor ion ve- 
locities are less well measured and only contribute a small 
fraction, so using the proton velocity for v is a reasonable ap- 
proximation. 

Example spectra of E^, E\, and Ei-, calculated over a 2 day 
intei-val (6th July 2007 00:00:00 - 8th July 2007 00:00:00) 
using the kinetic normalization are shown in Figure [T] Over 
most of the frequency range E\, is higher amplitude and 
steeper than Ey. Above « 3 x 10"^ Hz, Ey becomes larger 
than E\, due to noise in the velocity data (see Appendix). All 
spectral quantities were measured over the range 10~^-10~^ 
Hz, which is at low enough frequencies to avoid velocity noise 
affecting the analysis (see Appendix), at high enough frequen- 
cies to avoid the outer scale of the turbulence and gives a 
large enough range to obtain reliable spectral indices. This 
range corresponds to typical length scales of ^ 5 x lO'* km to 
'-^ 5 X 10^ km under the iTavlorl (11938b hypothesis, which is 
smaller than the typical correlation length and larger than the 
typical proton gyroradius. 

It can be seen in Figure [T]that the residual energy spectrum 
is noisy, in particular above 10"^ Hz but also below, since it 
is the difference between two similar spectra. Several meth- 
ods were attempted to measure its scaling: (a) a linear least 
squares fit in log space of E^, (b) a linear least squares fit in 
log space of the difference between the linear least squares 
fits in log space of Ey and E\,, (c) a linear least squares fit in 
log space of Er calculated using a wavelet transforms of v and 
b. Methods (a) and (c) were found to produce a large scat- 
ter in residual energy spectral indices since they are fits to E^, 
which is noisy. While Ey, E\, and E^ cannot, in general, all be 
power laws, in practice they approximately are, so method (b) 



CHEN ET AL. 




-0.5 



0.5 1 



Figure 2. Scatter plot of normalized residual energy ay and normalized 
cross helicity CTc using (a) MHD normalization and (b) kinetic normalization. 
Mathematically possible values lie within the circle. 



Table 1 

Mean values in the MHD and 
kinetic normalizations 



Parameter MHD Kinetic 



rA 



-0.43 
0.40 
0.40 

3.37 



-0.19 
0.46 
0.71 
4.45 



was found to give the most accurate scaling of Ei. In the rest 
of this analysis we use the results from method (b). Fits of 
Ey and E^, are marked on Figure [1] along with the difference 
of these fits, which can be seen to be almost power law and 
match £, well. All other spectral indices were calculated with 
linear least squares fits. 

In addition to the spectral indices of Ey, E^, E+, E-, E^, E^, 
El, and the average values of CTi, CTc, ta and r^, three other pa- 
rameters were calculated in each interval so that variation with 
solar wind conditions could be studied. These were the mean 
solar wind speed, Vsw, the normalized fluctuation amplitude, 
5B/B = a/(|(5B|2)/ ^|g|^^ averaged over the fit range, and the 
proton collisiqnal ag e, A^, obtained from the SWE parameters 
( lMarucaet"ani2011l) . 

3. RESULTS 
3.1. MHD vs Kinetic Normalization 

As described in Section |2] the analysis was performed with 
both the standard MHD normalization of the magnetic field 
(Equation [U and the kinetic normalization with temperature 
anisotropics and species drifts (Equation |2). Scatter plots of 
the normalized residual energy, (jj, and normalized cross he- 
licity, (Tc, (imbalance) for both normalizations are shown in 
Figure|2] The definitions of these quantities require a^ + al < 
1 (.Bavassano et al._1998.) . i.e., the points are constrained to lie 
within the dashed circles. 

It can be seen that the majority of points (« 80% in both 
cases) have positive values of CTc, corresponding to domi- 
nant propagation away from the Sun, as is well known (e.g., 
|Belcher& Davis 1971). Periods of Cc < may arise natu- 
rally from the turbulence or could be due to incorrect recti- 
fication due to polarity inversions (iBalogh et al.lll999b . The 
mean values of CTc are 0.40 with the MHD normaUzation and 
0.46 with the kinetic normalization. A more significant dif- 
ference can be seen in the a^ distribution, a^ is systemati- 
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Figure 3. PDFs of Alfven ratio r^ and Elsasser ratio r^ using the MHD 
normalization (dashed lines) and kinetic normalization (solid Hnes). 

cally more negative in the MHD normalization compared to 
the kinetic normalization, having mean values of -0.43 and 
-0.19 respectively. While the residual energy is still, on aver- 
age, negative, as predicted by turbulence theo ry (Wan g et aTl 
1201 U lBoldvrevetal.l2012bl: iGogoberidze et al..2012ah using 
the kinetic normalization more appropriate to the solar wind 
results in significantly less residual energy. 

The more natural variables in which to examine this dif- 
ference, perhaps, are the Alfven and Elsasser ratios, directly 
related to CTi and CTc through Equations |6]-|9] From Figure |3] 
it can be seen that r^ and te are approximately log-normally 
distributed and there is a much broader distribution of te than 
ta. Both are more positive with the kinetic normalization, 
and their mean values, calculated in log space, are r^ = 0.40 
and te = 3.37 in the MHD normalization and r^ = 0.71 and 
te = 4.45 in the kinetic normalization. These values are sum- 
marized in Table[T]and the standard error of the mean for all of 
these quantities is between 0.5% and 2%. This shows that r^ 
in the solar wind is significantly closer to unity than previous 
measurements suggested, but it is still < 1 at 1 AU. 

It should be pointed out that a small part of the differ- 
ence between the MHD and kinetic normalizations is due 
to the different instruments used to measure particle density 
(3DP/PESA-L for the MHD normalization and SWE/FC for 
the kinetic normalization), with a measured, rather than es- 
timated, alpha particle density in the kinetic normalization. 
This was done so that the MHD normalization can be com- 
pared to previous measurements. Calculating the results using 
the SWE/FC densities in the MHD normalization, however, 
shows that the majority of the difference is due to the tem- 
perature anisotropy and drift corrections. For the rest of this 
section, the kinetic normalization is used, since we consider 
the SWE/FC densities to be more accurate and the kinetic nor- 
malization more appropriate to the solar wind. 



3.2. Spectral Index Distributions 

Histograms of the spectral indices of E^, E^, Ei and E^ 
are shown in Figure |4] The distributio ns of E\,, Ey and 
El indices are similar to previous studies (ISmith et al .112006 



Vasauez et all 12007 
Boldvrevetal.11201 



I 



Tessein et al.ll2009t IChen et all 1201 la : 



Borovskvll2012h . with mean values 



1.69, -1.52 and -1.62 respectively. Due to the large number 
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Figure 4. Histograms of tlie measured spectral indices in the spacecraft fre- 
quency range 10"'-10"- Hz. The black dotted lines con'espond to various 
spectral index predictions. 

of intervals, the standard error of the mean is small, 0.002 
in each case. The spread of the £■, spectral indices is larger, 
which reflects the fact that this quantity is noisier since it is 
the difference between two similar spectra. It can still be seen, 
however, that the mean value of -1.91 is significantly steeper 
than the other three spectr a, as predicted by some resid ual en- 
ergy th eories (iMiiUer & G rappin 2004, 2005; Boldyre vetaTI 
l2012al) . The E+, E- and E^ spectra (not shown) are similar 
to that of the total energy, having mean values of -1.62, - 
1.61 and -1.60 respectively. The linear correlation coefficient 
between the E\, and E^ spectral indices is 0.44 ± 0.03 and be- 
tween the E^ and E^ spectral indices is 0.03 ± 0.03, where the 
uncertainty is the 99% confidence interval. This suggests that 
the properties of the residual energy are related to the mag- 
netic field, rather than the velocity. 

3.3. Dependence on Parameters 

The large data set enables the variation of the spectral in- 
dices with solar wind parameters to be studied. Figure |5] 
shows the mean spectral indices as a function of solar wind 
speed v'sw, magnitude of imbalance |crc|, collisional age Ac and 
fluctuation amplitude 5B/B. The error bars are ± 2 standard 
errors of the mean. It can be seen that the velocity spectral 
index remains close to -3/2 for for all values of these param- 
eters, but there are systematic trends in the other three. The 
magnetic field and total energy follow each other as expected, 
since the magnetic energy generally dominates. The residual 
energy also generally follows the trends of E^ and Ei but is 
overall steeper The spectral indices of E+, E- and E^ are the 
same to within errors as the total energy E^ in nearly cases, so 
are not shown. The only exception is that E- becomes shal- 
lower than the other spectra for large \ac\, large Vsw and small 
Ac. This is an artifact due to E- being low amplitude and close 
to the noise floor at large imbalance, which occurs in fast, low 
collisionality wind. 

The variation of spectral indices with Vjw and |(Tc| appears 
similar, which is likely due to the underlying positive cor- 
relation between solar wind speed and imbalance. Exam- 
ining the Eb spectral index as a function of both Vsw and 
IctcI indicates that a ^ is th e main cause of this variation. 
IPodesta & Boroys ky, ( 120101) also concluded that CTc is the pa- 



rameter controlling the variation of total energy spectrum. A 
similar explanation can be used to describe the Ac depen- 
dence, since Ac oc v^^. Again, examining the E^ spectral index 
as a function of Ac and |crc| shows Cc to be the main cause of 
this variation. It makes sense that a^, rather than Ac, con- 
trols the spectrum, since Alfv enic turbulence is not thou ght to 
depend on the collisionality (ISchekochihin et al]l200^ . The 
question then becomes what is causing the JctcI trend. At 
I CTc I = 1, the residual energy is constrained to be CTi- = so 
Ey and E\, take the same scaling, which can be seen in Figure 
ISj). As jcTcl decreases, the residual energy becomes non-zero 
(Tr < and Eb becomes steeper than E^, which is consistent 
with the properties of the residual energy (iMuller & Graopijil 
l2004l2005l:lBoldvrev et al.l2012aH Gogoberidze et al. 2012ai). 

The residual energy has a spectral index between -2 and - 
1.9 for < JCTcl < 0.75 (Figure [Sj)). This is consistent with 
the predicted value of -2 (Boldyrev et al. 2012a) and recent 
simulations of strong MHD turbulence, which cov er a similar 
range of imbalance and have a -1.9 spectrum (Bo ldyrev et al.l 
1201 lb- This scaling is not consistent with the theory of 
Miiller & Grappin (2004, 2005), for which the prediction is - 
7/3 (since the total energy spectrum is -5/3 over this range of 
|ctc|), or the theory of Go goberidze et al. (2012a), for which 
the prediction is -5/3. For |(Tc| > 0.75 the residual energy 
spectrum becomes shallower One possibility for this is that 
the highly imbalanced turbulence is less evolved due to its 
weaker nonlinear interaction and the residual energy has not 
had time to develop its steady state spectrum. A more mun- 
dane, but perhaps more likely, explanation is due to inaccura- 
cies in the normalization. At large JctcI the residual energy is 
small and the E^ and E^ spectra should be of similar ampli- 
tude, so a small error in the magnetic field normalization will 
cause a large error in the residual energy spectrum, biased to- 
wards shallower spectra. 

The variation of spectral indices with 5B/B (Figure |5}l) 
shows a slight trend towards shallower spectra at low am- 
plitudes. Examining the E\, spectral index as a function of 
both 6B/B and JctcI shows that both variables separately cor- 
relate to the spectral index. While the amplitudes are lower 
at small SB/B, the magnetic field spectra are well above the 
magnetometer noise floor, so the trend in Ey, appears to b e 
physical. This is consistent with the result of iLi et al.l (120111) . 
in that lower amplitude turbulence would be expected to pro- 
duce fewer strong discontinuities. 

3.4. Local Properties 

Figure |6k-d sh ows the continuous wavelet spectrograms 
(iTorrence & Com po 1998) of normalized cross helicity CTc and 
normalized residual energy CTr for two of the intervals. Such 
spectrograms have been used previously to stud y the nature of 
large scale fluctuations (iLucek & Baloghll998i) and flux ropes 
(Telloni et al. 2012) in the solar wind. The two intervals were 
chosen to be globally balanced (ctc = -0.13, te = 0.81) and 
globally imbalanced (a^ = 0.98, te = 95), covering the range 
of typical intervals. 

It can be seen in Figure |6^ that the globally balanced 
interval is made up locally imbalanced patches at all scales as 
expected for Alfvenic turbulence (£.g., Matthaeus et alj 
2008: iPerez & BoldvrevI l2009t [Boldyrev et al. 20091 
iPerez & BoldyrevI I2010I) . Such self-organization could 
be due to the nonlinear interaction in imbalanced patches 
being reduced so that they pers ist longer (|Per ez et al. 20121 
or to a rapid alignment process (iMatthaeus et al.ii2008i) . This 
interval also has local patches of both positive and negative 
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Figure 5. Dependence of spectral indices on (a) solar wind speed Vsw , (b) magnitude of imbalance | CTc | , (c) coUisional age Ac and (d) fluctuation amplitude 5B/B. 
The eiTor bars represent 2 standard errors of the mean. The black dotted lines coiTespond to different spectral index predictions. 



residual energy (Figure |6];), which is consistent with MHD 
turbulence, since the prediction of negative residual energy 
is for a statistical ensemble only, not individual fluctuations 
(IWang et aL.201 1). This figure shows that solar wind turbu- 
lence does contain localized patches where both velocity and 
magnetic energy dominate, while having a globally negative 
residual energy (ci = -0.20, ta = 0.70). This can also be 
seen in MHD turbulence simulations (^P erez & BoldyrevI 
120091 Figure 2), in which there are regions where either 
the magnetic field or velocity dominates. As expected, the 
imbalanced interval consists of large CTc values (Figure |6j5), 
which limits the spread of a^ values (Figure |6}i), due to the 
relation a'^ + al < 1 (Bavassan o et al.lll998l) . 

It is important to determine whether the local patches in 
Figure |6] are physical or due to instrumental noise. It can be 
seen from Figure[T]that the velocity noise is ^ 60 km^s"^Hz"' , 
whereas the total energy spectrum, over the range 10~^-10~^ 
Hz, varies from ^ 3 x lO"*- 1 x 10^ km^s"^Hz"\ which is a 
factor of ^ 40-2000 larger than the noise. The fluctuations in 
Figure|6]are on average around half the total energy, so cannot 
be due to this noise. 

Also shown in Figure |6^-j for each interval are the PDFs 
of ac, (Jr, fE, rp, and alignment angles cos(6'+_) = 5z^ ■ 
5z-/(|5z+||(5z-|) and cos(6i,,fc) = (5v • (5b/(|(5v||(5b|), calculated 
from 2-point differences: (5x = x(f)-x(f + t), where l/r = 



3 X 10 Hz. The a^ and Ci distributions (Figure |6^,f) con- 
firm the findings from the wavelet spectrograms, that the bal- 
anced interval is made up of patches of local imbalance and 
non-zero residual energy, and the imbalanced interval has a 
strongly peaked Cc distribution and narrower a^ distribution. 
The PDFs of rp, and r^ (Figure |6^,h) are approximately log- 
normally distributed, with standard deviations of 0.63 for r^ 
and 0.65 for ta in the balanced interval and 0.60 for te and 
0. 17 for ta in the imbalanced interval. This shows that imbal- 
anced turbulence has a similar spread of patches of local im- 
balance as balanced turbulence, but the mean value is larger, 

so the spread of residual energy patche s is smaller. 

For comparison to previous studies (Matthaeus et al."2008i 
jPpdes ta et al. 2009; Perez & Boldyrev 2009; Osman et af] 
120 llh . the PDFs of the cosine of the alignment angles are 
also shown (Figure |6j,j). The black dashed lines show the 
PDF for an even distribution of angles 6 and the black dot- 
ted lines show the PDF of the angle between two randomly 
oriented vectors. The behavior of cos(6'fei,), which is re- 
lated to the imbalan ce, is roug hly consistent with previous 
measurements (Matthaeus et al. 2008; Podestaetal. 20091 
lOsman et al.ll201 ll) 7having a broad distribution for small im- 
balance and becoming strongly peaked for large imbalance. 
cos(0+_), which is related to the residual energy, has a similar 
distribution to cos(6'/„.) in the balanced interval, showing that 
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Figure 6. Local properties of globally balanced (left panels) and globally imbalanced (right panels) turbulence, (a-d) Wavelet spectrograms of (Tc and (Tr. (e-j) 
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the turbulence also consists of a range of alignment angles of 
the Elsasser variable fluctuations. 

4. DISCUSSION 

In this paper, we have shown that when the magnetic 
field is converted to velocity units using a kinetic normal- 
ization that includes pressure anisotropies and drifts, the ra- 
tio of energy in velocity fluctuations to magnetic fluctua- 
tions is closer to unity, with a mean value of ta = 0.71, 
compared to ta = 0.40 with the standard MHD normaliza- 
tion. The residual energy, i.e., the difference in energy be- 
tween velocity and magnetic field fluctuations, is statisti- 
cally negative (although contains local patches of both signs 
at all scales) and has an average spectral index of -1.91. 
These properties are consistent with theoretical descriptions 
and numerical simulations of MHD turbulence ([Gr appin et al. 



T983';'Muller & Grappin'2004','2005\'Boldvrev &Perez 2009i 
Perez & Boldvrev 2009; Wang et al. 2011; Boldvrev et aQ 
201U ,2012ab; Gogoberi dze et al.ll2012aD . 

The magnetic field has an average spectral index 
close to -5/3, consistent with previous large studies 
dSmith et al.ll2006t IVasquez e t al.ll2007 "bt iTessein et al .ll2009l: 
Podesta & Borovskv 2010; Ch en et al.l l2011a; Boldvr evetaLl 
201 It iBorovskvi |2012|) . The velocity spectral index was 



found to be close to -3/2, as in previous studies with the 



Wind (Mangenev et al. 2001; Podesta et al. 2007; Sale niet alj 
[2009; .Podesta & Borovskv. .20 10; .Boldvrev et al...201 ll) and 
ARTEMIS JChenetaLl 1201 lab spacecraft. ACE measure- 
ments, howev er, produce a shallower spectrum, closer to -1.4 
( ITessein et al . 2009; Boldyrev et al. 2011; Borovskv 2012). 
The reason for this difference is not known, but could be 
related to the lower frequency range used in the studies of 
iBoldvrev et"an (120111) and Bor ovskv! (120 12h being affected by 
the outer scale and to high frequenc y noise (see Appendix) 
in the study of ITessein et al.l (120091) . The average residual 
energy spectral index of -1.91 found in this paper is steeper 
than the tw o values quoted in earlier studies (Tu et al. 19891 
iPodesta &B orovskv 2010). This is likely due to our use of 
the kine tic n ormalization more appropriate to the solar wind 
(Section lTTj i; with the MHD normalization, the mean residual 
energy spectral index with this data set is -1 .79. 

The residual energy in the solar wind fluctuations can be in- 
jected by the driving or can arise naturally from the turbulent 
interactions. An interesting area of future study would be to 
determine how much is due to the driving and how much is 
due to the turbulent dynamics, perhaps by a detailed compar- 
ison to simulations of MHD turbulence. 



This work was supported by NASA contract 
NNN06AA01C and NASA grant NNX09AE41G. Wind data 



CHEN ET AL. 



Podesta et all 12 0061: iPodesta & Borovskvl 120101: IWicks et alj 
201 U iGogoberi dze et al.ll2012bh . This flattening can be seen 
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Figure 7. (a) Power spectra of magnetic field E\, and velocity Ey from Figure 
[T], (b) quantization noise added to Eb, (c) white noise added to Ei,, (d) model 
-5/3 spectrum with and without aliasing. Spectral indices measured between 
10"^-10"~ Hz are given in the legends. 
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APPENDIX 

INSTRUMENTAL NOISE 

To measure accurate spectra, it is important to ensure that 
the fluctuations are not affected by instrumental noise that 
would cause a systematic error on the amplitudes and spec- 
tral indices. For the amplitudes considered in this paper, 
noise in the magnetic field measurements is insignificant. 
The noise on the velocity measurements, however, requires 
a more careful analysis. It is well known that various types 
of noise, here defined to be anything other than the physical 
values, can affect turbulence spectra by causing flattening at 
high frequencies (.Russell , 1972; Johnstone & KraukliS ;1998i: 



in the velocity spectrum in Figure [T] For the velocity mea- 
surements there are three major sources of such noise in the 
fluctuations (here we do not consider DC offsets): quantiza- 
tion noise due to limited telemetry (Russell 1972), Poisson 
noise due to the discrete nature of particles entering the de- 
tector dJohnstone & Krauklislll998[). and aliasing du e to dis- 
creet sampling (iRusselll 119721: IPodesta etall 120061) . Here, 
we consider each of these in turn, adding them to the mag- 
netic field s pectrum to inv estigate their effect, in the spirit of 
IGogoberidz e et al.l (I2012bi) . but considering physical sources. 

Quantization noise. Figure |7^ shows the E^ and Ey spectra 
from Figure [T] To determine the level of quantization in the 
velocity moments, a PDF of velocity differences at the data 
resolution (w 3 s) was calculated for each component. The 
PDF shows peaks at multiples of 2.3 km s~' for the GSE x 
component (the Sunward direction) and multiples of 1.3 km 
s"' for the other components. These were taken to be the 
quantization levels, although there is some variation about 
these values. Figure |7j) shows Ey as before along with E\, 
normalized to Ey at 10"^ Hz, calculated from b quantized to 
the same level as the velocity measurements. It can be seen 
that this quantization does not significantly affect the spec- 
trum and the spectral index remains unchanged in the range 
lO'-'-lO"^ Hz. This suggests that quantization noise is not the 
main source of the flattening of the velocity spectrum at high 
frequencies. 

Poisson noise. Each velocity moment is calculated from 
a measured distribution function, which is typically made up 
of a few thousand particle counts. This finite number leads 
to Poisson noise, which causes a white noise spectrum in the 
velocity moments, with an amplitude that is s pecific to the 
instrument design dJohnstone & Krauklislll998h . To test this 
effect, Gaussian noise of standard deviation 1.6 km s"' was 
added to the magnetic field spectrum to cause it to become 
equal to the velocity spectrum at high frequencies (Figure|7};). 
This represents the maximum amount of Poisson noise (in re- 
ality it may be smaller than or equal to this) and is within 
an order of magnitude of that found byl Johnstone & KrauklisI 
(11998). It can be seen that even with this maximum amount 
of Poisson noise, the spectral index remains relatively un- 
changed at -1 .64, compared to -1 .66. This suggests that while 
Poisson noise may account for the high frequency flattening, 
it is not the cause of the difference in spectral indices between 
Ey and E^. 

Aliasing. The effect of aliasing on spacecraft power spectra 
was discussed by Russell (1972) in the context of magnetic 
field measurements. A power spectrum P{f) sampled with 
a Nyquist frequency of /n will result in a measured power 
spectrum of 

R 

P'(f) = Pif) + Yl [P(2nfy-f) + P(2nfy + f)] , (Al) 

where 2R+1 is the ratio of the instrument upper cutoff fre- 
quency to the Nyquist frequency. The 3DP instrument on 
Wind generates an onboard velocity moment on each spin of 
the spacecraft, resulting in a Nyquist frequency of /n = 0.16 
Hz. Since the bulk velocity is larger than the thermal veloc- 
ity, the solar wind is only seen by the detector for a fraction 
of the spin period, which can be used to determine the upper 
cutoff frequency in this interval: the solar wind speed is 440 



km s"', the proton thermal speed is 39 km s"', so the 3 sigma 
solar wind beam width is in view for 6% of the spin, i.e., 0.18 
s, giving an upper cutoff frequency of 2.7 Hz. This gives a 
ratio of 16.7 and /? sa 8. Figure|7}l shows a -5/3 spectrum to- 
gether with the aliased version calculated from Equation lAll 
with these numbers. To construct this, it was assumed that the 
-5/3 spectrum continues up to 2.7 Hz; this is a worst case sce- 
nario, since the ion velocity spectrum has been measured to be 
steeper than this above 0.2 Hz (Safranko va et al.l[2013l) . It can 
be seen that the aliasing can cause some high frequency flat- 
tening but does not significantly alter the spectral index from 
10-3-10-2 Hz. 

We have shown that these three noise sources are individ- 
ually not likely to be affecting the velocity measurements in 
the range 10~^-10~2 Hz, but what about a combination? It 
is possible to make a model spectrum with the same spectral 
index as Ey in this frequency range by taking a -5/3 spec- 
trum, adding white noise at 70 km^s-^Hz-' then calculating 
the aliased spectrum. This, however, produces far more flat- 
tening at high frequencies than is seen in the Ey data (flat- 
tening to > 10^ km^s-^Hz-') and results in a non-power law 
spectrum in lO-^-lO-^ Hz, which is not observed. 

To conclude this Appendix, Poisson noise and aliasing are 
likely to be causing the velocity spectrum to become artifi- 
cially flat at high frequencies, but not significantly affecting 
the measurements in the range lO-^-lO-^ Hz that is consid- 
ered in this paper There are a host of other so urces of er- 
rors for particle instruments (IWuest et al.l 120071) . which are 
expected to be smaller than those discussed here for the cur- 
rent measurements, but could be considered in a more detailed 
analysis. 
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